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1. Introduction 

Consider an elastic medium which occupies a bounded domain Q C M 3 with smooth 
boundary dfl and exterior normal v. Displacement is a time-dependent vector field 
u(t, ■) on Q. Small displacements satisfy, in a source-free medium, the equations for 
(linearized) elastodynamics, 

(1) pd 2 u/dt 2 = V ■ S with S = R + VuR + CE. 

Here < p G C°°(Q) denotes the density, S is the Piola-Kirchhoff stress tensor which 
obeys the relation SF T = FS T where F = I + Vw is the deformation gradient. 
Divergence and transpose are taken with respect to the Euclidean metric | • |. The 
elasticity tensor C maps infinitesimal strain tensors E = (Vw + Vw T )/2 to symmetric 
stress tensors CE. C represents the material properties on the elastic medium. R(x), 
the residual stress tensor, is a symmetric 3 x 3-matrix, C°° on Q. It satisfies V-i? = 0. 
See ||Cur84 Sect. 23], ||Hog86|| . We call 



(2) Pu = -p d 2 u/dt 2 + V • (R + VuR + CE) 

the operator for elastodynamics. P is isotropic if the elasticity tensor is as follows, 

(3) CE = A tr(E)I + 2pE with < A, p e C°°(ty. 

A and p are the Lame parameters of the elastic medium. 

The inverse problem for operators of elastodynamics is to recover as much as pos- 
sible of the elasticity tensor C and of the residual stress tensor R from measurements 
performed at the space-time boundary K x dfl. See ||ML87|| for the beginnings of 



an acoustoelastic theory of residual stress determination based on wave propagation 
methods. 

We deal with an inverse problem for subclasses C(L,e) of operators for isotropic 
elastodynamics. Here L, e > 0, and, by definition, P e £(L,e) if and only if 

(4) X(x) + 2p(x),l/p(x),l/p(x) < L when x e H 
and 

(5) — when x E Q. 



The first author thanks the Department of Mathematics at the University of Washington for 
its hospitality during his visit in fall 2000. The second author is partly supported by NSF grant 
DMS-0070488 and a John Simon Guggenheim fellowship. The second author also thanks MSRI 
for partial support and for providing a very stimulating environment during the inverse problems 
program in fall 2001. 

1 



2 



SONKE HANSEN AND GUNTHER UHLMANN 



If e > is sufficiently small the initial boundary value problem is well-posed and 
microlocal parametrices exist. Assumptions like (||) with e small have been introduced 
before to ascertain well-posedness. See, e.g., |[Rob97|| for the static case. 
The (hyperbolic) Dirichlet-to-Neumann map 



(6) A : u\ Rxd n i-> v • S\ 



IxdQ 



encodes boundary measurements. Here u solves ([!]) with zero initial data. We say 
that a property of an operator for elastodynamics from a given class is determined 
by boundary measurements if the property is the same for any two operators in the 
class with identical Dirichlet-to-Neumann maps. 

A useful approach to inverse problems consists in using high-frequency waves u 
generated by boundary data with singularities. From travel times of singularities 
of u recorded at dQ one then aims to recover the requested properties. The latter 
problem is called an inverse kinematic problem. Obviously, an important step is 
to prove that travel times are in fact determined by boundary measurements. The 
main goal of this paper is provide a result of this kind which is applicable also when 
caustics may develop. 

We study the propagation of polarization in the sense of Dencker |Den82|| for the 



initial boundary problem of the operator for isotropic elastodynamics, P. In Propo- 
sition |0] we show that P is a system of real principal type if the residual stress R 
satisfies 

(7) Mx)|£| 2 + R{x)Z ■ f > when (x,£) G T*(H) \ 0. 
If © holds then 

(8) (U)s= (^m 2 + R(x)^^/p(x), 

(9) (£, Op = ((A(x) + 2Mx))|£| 2 + R(x)£ ■ £)/p(x), 

(x,£) G T*(f2), are the duals g~gj P of Riemannian metrics gs/p on Q. The character- 
istic variety of P is the union of the subvarieties r 2 — (£, £) g = and r 2 — (£, £) p = 
which correspond to shear and compressional waves, respectively. 

The lens map or scattering relation S of a metric g on Q is defined as follows. 
Consider bicharacteristic curves, 7 : [a, b] — * T*(Q x R), of the Hamilton function 
H(t, x, t, £) = r 2 — g^ 1 (x, £) which satisfy the following: jQa, b[) lies over the interior, 
7 intersects the boundary non-tangentially at 7(a) and 7(6), and time increases along 
7. By definition, S is the subset of (T*(R x dQ ) \ 0) 2 obtained by projecting end- 
point pairs (7(6), 7(a)). It is well-known that S is a homogeneous canonical relation 
on T*(R x dQ) \ 0. (See [ Gui77 ] for the concept of a scattering relation.) S is a 



diffeomorphism between open subsets of T*(R x dQ) \0. We denote by S$ (resp. Sp) 
the lens map of gs (resp. gp) and call it the shear (resp. compressional) lens map. 
Our main result is the following. 



Theorem 1.1. Given L > there exists e > such that in the class C(L,e) the 
shear and the compressional lens maps are determined by boundary measurements. 
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Note that travel times of shear and compressional waves are recovered separately 
from boundary measurements. 

Let g be Riemannian metric on Q. Denote by D the open subset of dfl x dQ 
which consists of the pairs (x, y) of boundary points which can be joined by a ge- 
odesic which passes through the interior except for the endpoints x and y where it 
intersects dfl transversally. By definition, the boundary distance function of (Q, g) 
is the function d : D — > [0, oof which assigns to (x,y) G D the geodesic distance, i.e., 
the infimum of the lengths of such geodesies. If (fl, g) is strictly convex then D is 
the complement of the diagonal and d is smooth. Geodesies of g are projections of 
bicharacteristic curves of t 2 — <7 -1 (a;,£) = 0. Geodesic distances equal travel times. 
When d is smooth it is a generating function of (a subset of) the lens maps S of g, 
i.e., ((ti, 3?i, t, ( t , x , t, g o)) G <S if h - t = d(x 1 ,x Q ) and ^ = -rdd(x 1 ,x )/dx j 
for j = 0,1. (See ||Gar35| 1 , ||GS77|| .) Clearly, S determines d. Hence we have the 



following corollary of Theorem Here we call the boundary distance functions of 
the metrics g$ and gp the shear and the compressional boundary distance functions 
ds and dp, respectively. 

Corollary 1.2. Given L > there exists e > such that in the class C(L,e) the 
shear and the compressional boundary distance functions are determined by boundary 
measurements. 

Rachele ||Rac00b|| has a similar result under additional assumptions which exclude 



conjugate points. Note that our result allows the presence of conjugate points. 

In the case R = the metrics g$ and gp are conformal to the Euclidean metric. 
Mukhometov [|Muk82|| solved the inverse kinematic problem for conformal classes 



of metrics under assumptions which exclude conjugate points. Corollary [L2| and, 
e.g., Croke's theorem [ Uro9 1 , Theorem C] imply the following uniqueness result of 
Rachele. 



Corollary 1.3. [RacOOa] , Theorem 1] In the class of operators of isotropic elasto- 



dynamics with vanishing residual stresses, and with (fl } gs), (Q,gp) strictly convex, 
the compressional speeds and shear speeds, cp = a/(A + 2/i)/p and cs = y/jjjp, are 
determined by boundary measurements. 

If residual stresses do not vanish the metrics become anisotropic. From Corol- 
lary [L2| and a result of Stefanov-Uhlmann on the anisotropic inverse kinematic prob- 
lem ||SU98| , Theorem 1.1] we deduce the following result. 

Corollary 1.4. There is a C 12 (Q) neighbourhood U of the euclidean metric such 
that the following holds. Let and be operators of isotropic elastodynamics. 
Assume = . Assume fl strictly convex with respect to the metrics and 
gp\ If gg\g^P € U then gy = ^f* s gf\ g^p = ^pgp 2 ^ with diffeomorphisms ^ s ,^p '■ 
fl — ► Q which leave the boundary fixed, i.e., \&s(x) = ^p(x) = x if x G dfl. 

In | jSU98| , Theorem 1.1] an additional flatness assumption at the boundary of Q is 
made. This assumption is superfluous in view of ||LSU01| , Theorem 2.1]. 



We prove Theorem |1.1| in section [| The facts needed about propagation of sin- 
gularities and polarizations in non-glancing boundary problems for systems of real 
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principal type are proved in section ^| for first order systems. These are applied to 
second order systems and to elastodynamics in sections |3| and £|, respectively. In 
particular, section [| contains an analysis of the Dirichlet-to-Neumann map A and its 
pseudo-differential properties. 



2. Singularities of First Order Boundary Problems 



We summarize some facts from the microlocal theory of boundary problems 

Dengg , Gerard | |Oer85|| , Melrose ||Mel8l|| 



results are due to Dencker 
Tay75 |. 



The 
and Tay- 



lor 

Let Z an open subset of half-space R + x M. n equipped with coordinates x > 
and y = . . . , y n ). £ and r\ = (rjx, ... , i] n ) are the dual coordinates in cotangent 
space. Denote the boundary and the interior of Z by Y and Z°, respectively. T>'(Z) 
denotes the space of extendible distributions on Z° . Pseudo-differential operators of 
order at most m on Y and on Z acting along Y are written A(y, D y ) G \l/ m (Y) and 
B(x,y,D y ) G ^r(Z), respectively. S m = S m (Y x W l ) and S™ = S m {Z x W l ) are 
the corresponding symbol spaces. Elements of are called tangential pseudo- 

differential operators. Symbols are always assumed polyhomogeneous (classical). 
Pseudo-differential operators will always be chosen properly supported. We denote 
by r Y u = u\ Y the restriction of u G V'(Z), when defined. 

We consider u G T>'(Z) K such that Pu G C°°(Z) K where P is a K x K system of 
pseudo-differential operators which are differential with respect to x, 



(10) 



P 



with P j = PJx,y,D y )e*i(Z), P = ld 



K ■ 



3=0 



The boundary wavefront set WFf, (it), defined in |Mel81|l , is a closed subset of the 
compressed cotangent bundle T*(Z). If the Pj are differential operators then Pu G 
C°° is a non-characteristic boundary problem and hence u is normally regular in the 
sense of Melrose ||Mel81| , II.9]. Recall from ||Mel81|| or ||Hor85| , 18.3] the following 
properties of a normally regular distribution u. u G C oo ([0,e[,V(W a )) locally near 
Y. Au is normally regular if A is a tangential pseudo-differential operator. The 
boundary wavefront set WF b (u) C T*(Y) UT*(Z°) C T*(Z). (y,r)) G T*(Y)\WF b (u) 
if and only if Au G C°°(Z) for some operator A = A(x,y,D y 



non-characteristic at (0,y,r]). 
intersection of the sets 



The polarization set WF 



pol 



G ^T( Z ) which is 
u) is, by definition the 



Ma = {(x,£;w) G T\Z°) x C K ; a{A){x^)w = 0} 



where 4g$° runs over all 1 x K systems such that Au G H^ S \Z°). See | Den82| | and 



Ger85|| for the precise definition and for results on the propagation of polarization 
along Hamilton orbits. 

Let P as in (|I0"D with principal symbol p. Following Dencker | Den82 , Definition 3.1] 
we say that P is of real principal type if, microlocally near a given point, the char- 
acteristic variety is given by q = with a scalar symbol q of real principal type 
and if there exists a matrix- valued symbol, p such that pp = qld.K- If we assume P 
of real principal type then H = H q is a Hamilton field of the characteristic variety 
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V = g _1 (0) of P. A point (y, rj) G T*Y \ is called glancing for P if, with respect to 
the natural projection, its preimage in V DTpZ contains a point where Hx = 0, else 
(y, rf) is called non-glancing for P. Bicharacteristics intersect the boundary transver- 
sally at non-glancing points. 



We now specialize to first order systems, 



m 



1. Let G = G(x,y,D y ) G ^]{Z) 



be an N x N matrix of tangential pseudo-differential operators with homogeneous 
principal symbol g. We assume that D x ld N — G is of real principal type. We are 
interested in the singularities of normally regular solutions of 



(11) 



D x w-G(x,y,D y )w = mod C°°(Z) 



N 



Let r]^) G T*(Y) \ non-glancing for D x ld N -G. 

The following decoupling lemma is due to Taylor |[Tay75|| in the case of simple real 
characteristics and to Gerard |Ger85|1 in the case of real principal type systems. 

Lemma 2.1. In a conic neighbourhood Y of(Q,y( \r]( '), the algebraic and geometric 
multiplicities of the real eigenvalues of g(x,y,r)) are equal and constant. There are 
homogeneous real-valued fi±, . . . , \ij G S\ which enumerate, in Y, the distinct real 
eigenvalues of g. Let Nj denote the multiplicity of fij. There is an elliptic N x N 
matrix S G such that microlocally near (0, y^°\ rj^), 

(12) (D x ld N -G)S = S(D x ld N -H) mod* t -°°. 

H G is a block matrix with non-zero entries only on the diagonal, 

( fii(x,y,D y )ld Nl \ 



(13) 



H 



^j(x,y,D y )ld 



Nj 



E 4 



V 



E- J 



The imaginary parts of the eigenvalues of the principal symbols of E + ,E_ G ty] 
positive and negative, respectively. 



are 



(o) „(o)> 



Proof. The following constructions hold in some conic neighbourhood Y of (0, y^°',rj 
Y may become smaller as the proof proceeds. 

Since A = D x ldw — G is of real principal type its characterictic variety is V — 
g _1 (0) with a scalar real principal type symbol q. The non-glancing assumption 
implies dq/d^ ^ at points (0, y^ - 1 , ^, 77^°-*) G V. By the implicit function theorem, 
the real eigenvalues of g(x, y, rj) are smooth homogeneous functions fii(x, y, rj) < . . . < 
/ij(x,y,r}) in Y. We extend them as homogeneous real valued symbols fii, ... G 
S}(Z x R n ). 

Let /i be a real eigenvalue of g(°> = g(0,y^°\r]^). We show that the geometric 
multiplicity of \i equals its algebraic multiplicity, 

(14) ker ((/i - g^Y) = ker (fi - ^ (0) ) , Vr G N. 

Let a = £ — g denote the principal symbol A. By the non-glancing hypothesis 
<9/<9£ is transversal to the characteristic variety det a = at (y^°\ r]^). The intrinsic 
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characterisation of real principal type |pen82| , Prop. 3.2] shows that da/d^ = Id maps 
the kernel of a isomorphically onto the cokernel of a at £ = fi. Hence 

(15) ker(/i-3 (0) ) n im^ - ^ (0) ) = 0. 

Equation (|14D easily follows from ([15]). 

Let 71, • • • ,7j be non-intersecting closed positively oriented Jordan curves in the 
complex plane such that 7j encloses fj,j(0,y ( -°\r]^) but no other eigenvalue of g^°\ 
(To enclose means that the winding number is non-zero.) 

(A - g(x, y, ^VM) -1 7^— 

is the spectral projector onto the sum of generalized eigenspaces associated with the 
eigenvalues enclosed by of g{x,y,rf)/\r]\. Clearly ker (fj,j — g^) C im-Kj. By ( |i4|) 
equality holds at (0, y(°\ r/°)). By [pen82| , Prop. 3.2] the dimension of ker — g^) 
is constant. Also the rank of 7ij is constant in T. Hence 

(17) ker(yUj — g) — im-Kj in T. 

It follows from (|17|) that the geometric and the algebraic multiplicities of the real 
eigenvalues of g coincide everywhere in T. Therefore we can find an elliptic N x N 
matrix s(x,y,r)) G S® such that s~ 1 gs G Si has, in T, the block structure of the 
principal symbol of the operator H claimed in (|13"D. 

Choose S with principal symbol equal to s. We obtain (|T2| ) with the error class 
replaced by ty®, however. We use the uncoupling technique of ||Tay75|| to obtain 
K £ ty^ 1 such that the error is if we replace S by S(ld+K). After doing this, 
however, H will only satisfy a weaker form than (13) with [ij(x, y, D y ) Idjy,- is replaced 
by [ij(x,y,Dy)l<i]s[. +Mj with some Mj e By | |Cer85| , Lemme 2.1.] there exist 
elliptic Nj x Nj matrices Ej e ty® such that 

((D x - fij(x, y, D y )) ld N . -Mj)Ej = Ej{(D x - /j,j(x, y, D y )) ld Nj ) 

holds modulo operators in ^>~°° . Let E G ty®, N x N, denote the diagonal block 
matrix with blocks E\,... , Ej and, in the lower right corner, Id. Finally, to remove 
the Mf's, we replace S by SE. □ 

Let B be a. K x N matrix in with homogeneous principal symbol b. Given 
h G V'(Y) K we wish to solve equation ( |TT| ) under the boundary condition specified 
by B and h, 

D x w -Gw = mod C 00 ^)^, 

(18) 

Sw| y = /i mod C°°(y) x 

Let M + UM_ = {1, . . . , J} be a disjoint union decomposing the set of real eigenvalues 
of g(0, y(°\ 77^) into two parts. We call the eigenvalue [ij forward (resp. backward) if 
j G M + (resp. j G M_). Correspondingly, we call characteristics and bicharacteristic 
curves forward or backward. In case D x Idjv —G is hyperbolic with respect to a time 
variable t(x, y) such a decomposition arises as follows. A bicharacteristic 7 issuing 
from the boundary into the interior is forward (resp. backward) if t increases (resp. 
decreases) along 7 
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We shall find a microlocal parametrix of the boundary problem ( |18|) if a condition 
of Lopatinski type holds. Define, for (y, 77) sufficiently close to (y^°\ 17^), the forward 
Lopatinski space as the following linear subspace of C^, 

(19) L+(y, 77) = im jf + (A - <?(0, y, r/))- 1 dA. 

7 + is a closed positively oriented Jordan curve in the complex plane which encloses 
the eigenvalues of = g(0,y^°\r]^) which are real and forward or which have 
positive imaginary part. 7 + encloses no other eigenvalues of g^°\ 

Proposition 2.2. Assume that b(0,y^°\rj^) maps L+(y(°> ,r)(°>) onto C K . Then 
there exists a conic neighbourhood T C T*(Y) of (y^- ' ,^ ') and an operator W : 
D'(Y) K — > V(Z) N such that the following holds. For every h G T>'(Y) K with 
WF(h) C T the distribution w = Wh G V(Z) is normally regular and solves 



(|18l). WF(w\z°) is contained in the union of the forward bicharacteristics which is- 
sue from WF(/i). Wo := ry o W G ty°(Y) is a N x K pseudo-differential operator. 
The principal symbol Wo(y,r]) of Wo maps C K into L^(y,r]) and satisfies bwo = Id in 

r. 

Proof. Let 5* and H as in Lemma |2.1|, and denote their principal symbols by s and 
h, respectively. Clearly, gs = sh. Hence L+ = sL^. The block structure of H and the 
partitioning into forward and backward eigenvalues defines a projector II on C N . U 
projects onto the subspace corresponding to the blocks fij(x, y, D y ) Id^., j G M + , and 
E + of H along the subspace corresponding to the blocks fij(x,y, D y ) Idjv,, j G M_, 
and E_ of H. Notice that L\ = HC . By assumption 

(20) C K = bL+ = bsL+ = bsUC N . 

The Cauchy problems 

(D x - N {x, y, D y ))v G C°°(Z), v\ Y G C°°(Y), 



are solved using scalar Fourier integral operators Vj, [ [Uui73| ]. The wavefront set of 
the solution v — Vjf is contained in the image of the bicharacteristics associated with 
£ — fij(x,y,rj) = which issue from WF(/). The parabolic system 

D x v — E + v G C°°(Z), v\y G C°°(Y), 



is solved using a Poisson operator V + , |[Lay75|| . The solution v = V + f has no sin- 



gularities in Z°. Therefore we may construct an operator V : T>'(Y) N — > T>'(Z) N 
such that the following holds for any / G T>'(Y) N . v = V f is normally regular, 
D x v-Hv G C°°(Z) , and WF(v\ z °) is contained in the union of the forward bichar- 
acteristics which issue from WF(/). Furthermore, modulo C°°(Y) N , v\y = II/. 

ryBSV is a K x iV system of pseudo-differential operators on Y . Its principal 
symbol bsll is, close to (y^°\ f/ ^), surjective by (p0|). Choose a N x K operator 
C G which is a right inverse, ryBSVC = Br Y SVC = Id. W = SVC satisfies 

the claims. □ 
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Remark 1. If the boundary data f is a Lagrangian distribution then the solution w = 
Wf is Lagrangian with respect to the forward characteristics. Rohrig |[Roh]| derives 



the transport equations for the principal symbol of w along the bicharacteristics. 

To prepare waves with specified polarization we need the following result about 
propagation of polarization at the boundary. Essentially this is a corollary of | Gcr85 . 
Theoreme 6.1]. 



Proposition 2.3. Letw £ T>'(Z) normally regular such that (y^°\ r/^) WFb(D x w- 
Gw). Assume w\y £ H^^^Y), s > 1. Let fi £ S 1 be a real eigenvalue of g(0, •) 
in a conic neighbourhood of (y^°\r)^). Let Q £ ty°(Y) with principal symbol equal 
to, in a neighbourhood of (y^ \ r)^ ') , the spectral projector on the eigenspace of the 
eigenvalue fi. Then (y(°\ rj^) £ WF^ s \Qw\y) if and only ifWF^ } (w) contains a 
Hamilton orbit above the fi-bicharacteristic which issues from (y^°\r]^). 



Proof. Choose a parametrix S^ 1 of S in Lemma and put w' = S -1 ^!. The 
hypotheses of the Proposition still hold with w replaced by w' and with G replaced 
by H of fll3|) . Let denote the projection to the components of the block which 
corresponds to fi in the block decomposition ([Hf). Then Q — £ ^^(Y) and, 
using the assumption on w\y, WF^(<5 M w|y) = WF^ s \Qw\y)- v = Q^w solves the 
diagonal system (y(°\ rj^) WFb(D x v — fi(x, y, D y )v ). The assertion follows from 
well-known results on propagation of singularities in the Cauchy problem for scalar 



strictly hyperbolic equations and from ||Den82| , Theorem 4.2]. □ 



3. Second Order Boundary Problems 

Here we reduce the Dirichlet problem for second order real principal systems to a 
boundary problem for a first order real principal type system. 

Let P = D 2 X Y&k +Pi(x, y, D y )D x + P2(x, y, D y ) be a K x K matrix of differential 
operators of second order. We are interested in the Dirichlet problem 

Pu = mod C°°(Z) K , 

(21) 

u\ Y = f mod C°°(Y) K . 



Any solution u is normally regular. 

We associate with (|2lf) an equivalent first order boundary problem as follows. 
Set N = 2K,G e J the N x N matrix 



(22) G 



(D y ) ld K 

-PliDy)- 1 ~P 1 



and B £ the K xN matrix with Bw — Wi, w — (wi, w^). Here (D y ) £ ty\ denotes 

1 /2 

the operator with full symbol (rj) — (1 + \r)\ ) £ Sj. 

Lemma 3.1. Let f £ V(Y) K and h = (D y )f. Solutions u of (|2T| ) andw = (w\,W2) 
of (0) are related as follows. If u solves (PI|) then w = (wi,W2) = ((D y )u,D x u) 
solves flTE|). Conversely, if w solves (|T^) then u = (D^^wi solves fl2"T|). 
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Proof. The first statement follows immediately from the definition of G and B. 
For the proof of the converse statement let w be a solution of fllSp. The first row 
of D x w = Gw and the ellipticity of (D y ) imply D x u = W2- Hence the second row 
implies Pu = 0. By our choice of B and h the boundary conditions are equivalent: 
(D y )u = Bw = h= (D y )f. □ 

Let (x, y, £, r/) G T*(Z), r\ ^ 0. Let p = £ 2 + pi£ + P2 denote the principal symbol 
of P. Then the principal symbol of G is 



9 

Lemma 3.2. Let f] ^ 0. Then 



\r)\ 
-p 2 /\v\ -Pi 



(23) (£-g')(£-g)= ( P p ) uAere </ = ( 



P \ J -I -pi -\v\ 

P2/\V\ 



and 



(24) kei ^-g)=^ jker(p). 

The characteristic varieties of P and D x ld.N —G are equal. If P is of real principal 
type then so is D x ld.N —G. 

Proof. Equation fl2"3] ) is verified by direct computation. Clearly, (£ — g)w = 
with w = (wi,W2), holds if and only if £w± = \rj\w2 and pw2 = 0. To prove the 
last assertion assume there is a K x K matrix of symbols, p, such that pp = gldx 
holds with a scalar real principal type smbol q. Then, using (p3|), we obtain a, N x N 
matrix of symbols, a, such that a(£ — g) = qldjy. □ 

Remark 2. Assume P of real principal type. LetCou = ((D y )u\Y, D x u\y) denote the 
Cauchy data of a solution of Pu = 0. It follows from Proposition \2. j| and Lemma \3.2j 
that WFp S j~ 1 ' 1 (u) contains a Hamilton orbit above a given bicharacteristic issuing from 

7 = (y, rj) G T*(Y)\0 if and only ifq G WF^ s \QC u) where Q G ty°(Y) with principal 
symbol equal to the spectral projector onto the eigenspace {((r))a,£a); p(0,y,C,,r])a = 
0} which corresponds to the given characteristic. 

We give sufficient conditions for the existence of a microlocal parametrix for the 
boundary problem (|21"|). 

Proposition 3.3. Assume P of real principal type. Let (y^ ' ,r)(°') G T*Y \0 be non- 
glancing for P. Let 7 + be a closed positively oriented Jordan curve which does not 
meet the poles of A t— > p(0, y(°\ A, r/^) and which has winding number 1 (resp. 0j 
respect to the poles with positive (resp. negative) imaginary part. Assume that 

(25) if > rank/ (A - g(0, w (0) , ?/ 0) )) 1 dX, 

(26) if < rank/ p(0, w (0) , A, j ? (0) )" 1 dX. 

J-y+ 
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Then there exists a conic neighbourhood T C T*(Y) of (y^°\ r/^) and an operator U : 
V'{Y) K -> V'{Z) K such that for any f E V'{Y) K with WF(/) C Y the distribution 
u = Uf E T>'(Z) K is normally regular and solves (pl|). WF(u\z°) is contained in 
the union of the forward bicharacteristics which issue from WF(/). U' := rYD x U E 
^/ 1 (Y) is a K x K pseudo- differential operator with principal symbol u' which satisfies 

(27) «V°W 0) ) I p(0,y^,X,r ] ^y 1 d\= [ Ap(0, y (0) , X, r]^)' 1 dX. 



7 t 



Proof. We use the equivalence, stated in Lemma |3.1| , of with the first order 
boundary problem (|I8|). 

A real eigenvalue of g(0, y(°\ r/°)) is, by definition, forward if it is enclosed by 7 + . 
First we show that our assumptions imply the following formula for the Lopatinski 

space, 

(28) L+(/>),r/°))=im( \ 1 \ lP } % ' 




From (|23|) we infer that the resolvent of g is 

(29) (A - gY 1 = ( * ^ P ^-i ) where p(A) = X 2 + Pl X + p 2 , X E 



* 



Ap(A)" 



Hence the right hand side in (|28f ) is contained in the left hand side. Equality follows 
from the dimension assumptions (p5|) and (EH). 



The principal symbol of B is 6 = (Id^-,0). Therefore (|28 ), (|25|) , and (f2q) imply 
= C x at (0, i/ -*, rp )). Proposition applies to give a solution operator of 



flTJ), W = (W U W 2 ) with Wi = BW. Define C/ = {D V )~ 1 W 1 {D V ) . It follows from 



Lemma g3] that PU = 04s, and r y C/ = Id, and W 2 (D y ) = D X U. Hence ((D y ),U') = 
ryW(D y ) E ^(Y). The principal symbol of ryW, (Id^, |77|~ «'(?/, rj)), maps into 
the Lopatinski space L^(y, if). Now we can read the formula ( p7|) off the equation 
The bound on WF(m|^°) follows from the bound on WF (w\z°) in Proposition 
□ 

4. Isotropic Elastodynamic Equations 

In following P denotes an operator for isotropic elastodynamics introduced in (0) 
and (|3|) such that (|7|) holds. 

The boundary problem Pu = g in R x f2, and m = on R x <9f2, has the variational 

o 

formulation: (pit, f) + a(w,t>) + (<7,i>) = 0, Vv E (H^(Q)) . Here a = ao + or, 
ao(M,f) = Jo tr (CE(u) E(v) T ) dx, and an(u,v) = f n tr ((Vw) -R (Vt> ) T ) drr. It fol- 
lows from Korn's inequality that a satisfies a coerciveness estimate |a (it, w)| > c||«||^ 
with a positive constant depending only on Q and on a lower bound on the Lame 
coefficient p, |DL76|| . Given L > there exists e > such that a is coercive if 
P E C(L,e). In fact, a# is absorbed into the coerciveness estimate if (|5]) is assumed 
with < e = e(L) sufficiently small. We use |pL76|, Thm. III. 4.1.] to conclude 



that the initial boundary value problem for P with Dirichlet boundary conditions in 



ELASTODYNAMICS WITH RESIDUAL STRESS 



11 



Hc S \M. x dfl) , s > 3, is well-posed. In particular, the Dirichlet-to-Neumann (DN) 
map (||) is defined, 

(30) A : H^ S+1] (R x dtt) 3 -> H {s) (R x <9ft) 3 if s > 2. 

Let (t, x,t, £) denote a generic point in T*(R x fi) \ 0. The Euclidean metric 
£ 2 = £ • £ is used to identify tangent and cotangent vectors of fl For t;, ( 6 C 3 the 
dot product is the analytic (non-Hermitian) extension, 77 • ( = rjid + 112(2 + V3(3- 
vr = 7r(£) = (£ <8> £)/(£' denotes the orthogonal projection onto a nonzero direction 
if £ G R 3 . 

As a consequence of (jfy the metrics defined in (H) and (|9|) satisfy 

(31) 0<(Z,Os<&Op if ^°- 



The norms associated with these metrics are denoted \£\ S / P = y (£,0s/p- 

Proposition 4.1. The scalar symbols qs/p(t,x,r,£) = p(x) (r 2 — (£, Os/p) an d th 
product qsqp are of real principal type. P is a system of real principal type with 
principal symbol 

(32) p = q s (Jd 3 -it) + q P ix. 

Proof. A straightforward computation gives the principal symbol p of P at (t, x,t,£) G 
T*(R x U) as follows: 

p = pr 2 ld 3 -(A + ® - ^ 2 W3 -(£ • Ida • 

Hence p has the asserted form. It follows from ( |3~i~D that and q P are of real principal 
type. Furthermore 

(33) q P {t,x,T,0<qs{t,x,r,0 if 0. 

Hence also g = gs^p is a scalar symbol of real principal type. Now q ,_1 (0) = 
(detp) X (0) and pp = qld for p = q P (ld 3 —n) + g 5 7T. According to |Pen82| , Def- 
inition 3.1] P is a system of real principal type with characteristic variety q = 0. 
□ 



Remark 3. Man ||Man98|| proposes for elastodynamics with residual stress R a more 
general constituitive law S = R + Vw • R + C-E 1 where the elasticity tensor C also 
depends linearly on R. In the isotropic case CE consists of the right-hand side in @ 
plus the R dependent terms 

(34) ft tr(£) tr(i2)J + (3 2 tv(R)E + ft ( tr(E)R + tr(ER)l) + ft + . 

In the inverse problem for real media the additional terms should not be neglected 
since typically R is much larger than the stress S. A straightforward calculation 
shows that the elastodynamic operator P with this isotropic stress-strain relation is 
still of real principal type in case ft = ft = 0, A + 2/i + ft tr(R) > p + ft ti(R)/2, 
and (p(x) + ft (x) tr(R(x))/2) |£| 2 + R(x, f ) > when (x, f ) G T*(H) \ 0. 
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We recall some notions of the microlocal theory of boundary problems and apply 
them to the system of elastodynamics. Let 7 = (t, x, r, £1) G T*(M x dfl) \ 0. This 
means that there is given (t,x, r, £) G T*(R x fi) with x G <9f2 and £1 = £|T x (an)- 7 
is called an elliptic, a hyperbolic, or a glancing point of S/P mode if the following 
quadratic equation in z, 

qs/p(t,x,T,£ - zv{x)) = 0, 

has no real roots, two distinct real roots, or a double real root, respectively. T*(K x 
dQ) \ decomposes into the disjoint union of the elliptic region £s/p, the hyperbolic 
region 7is/p, and the glancing hypersurface Qs/p of the S/P mode. Because of fl33|) 
we have £$ C £p and 7i P C H5. T*(R x dVL) \0 is the disjoint union of the hyperbolic 
region Tip, the mixed region Sp fl the elliptic region £3, and the glancing set 
Q = Gs U The lens maps satisfy Ss/p C 7Ys/p x T^s/p- 

A simple real root 2 is called forward (resp. backward) if the bicharacteristic curve 
starting in £ — zv enters 1R x Q when time increases (resp. decreases). Charac- 
teristics and bicharacteristics are called forward or backward correspondingly. Ob- 
serve from Hamiltons equations that a characteristic £ — zv, z real, of qs/p is for- 
ward (resp. backward) if r(£ — zv,v) s , p is positive (resp. negative). We denote by 
zs/p = Zs/p(t,x,r,£,u) the forward real root z or the complex root z with pos- 
itive imaginary part of qs/p{t,x,r,^ — zv) = 0. We shall use the abbreviation 
£s/p = £ - z s / P v(x). 

Given 5 > we define 

(35) T s = {(t,x,T,£\) £f(lx dQ) \0; |r| > 6\£\\}. 

Here |£i| = |£| if £ • z/ = 0. A is pseudo-differential microlocally at nonglancing points 
in rY 

Proposition 4.2. Let L, S > 0. Assume (|4]). There exists < e = e(L,5) such that 
under the assumption ([5]) the following holds. Given 7 G \ ^ i/ie DA'' map A equals, 
in a microlocal neighbourhood of (7, 7), a first order pseudo-differential operator with 
principal symbol given as follows. 

(36) ct(A) : a i-> X(a ■ £> + /i(a • z/)£' + //(£' ■ z/)a + (#£' ■ v)a 
when a G C£', resp. a - £' = ; £' = £p ; resp. £' = £5. 

We the need the following fact about the characteristics of P. 

Lemma 4.3. Let L,5 > 0. Assume (H). There exists < e = e(L,8) such that the 
following is true if (H) holds. For 7 = (i, x, r, £|) G r<5 \ we have z$ 7^ Zp, £s 2 7^ 0, 
£p 2 7^ 0, and 

(37) ( S '(p/ 0. 



We prove Lemma |4.3j in section [7|. 



Proof of Proposition |T2|. Choose e as in Lemma £D| Let 7 = (t, x, r, £|) G 



Ps \ (?. Flatten the boundary d£l near x with a change of coordinates such that the 
differential at x is orthogonal. 
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We consider the symbol p = p(t, x, •). Its inverse is p^ 1 = g^" 1 (Id3 — 7r) + q^n. Let 
7 + be a closed Jordan curve enclosing zs and zp but no other roots of q(r, £ — zv) = 0. 
Observe cs/p ■= (d/dz)qs/p(r,£ — zu)\ z=Zs/p ^ 0. By the residue theorem 

r d7 

(38) ' V V ' 2ttz 

= (4/c 5 ) (Id 3 -tt(&)) + (4/cp) vr(£p) 

for every non-negative integer j. We show that A is non-singular. Assume A w = 0. 
Then = £s • A w = (£s • £ p ) (£p • w)/{cpQ,) . Applying Lemma [O^ we infer £p • u> = 0. 



Therefore w = 7r(£s)n; and = (£p ■ £s)(£s • w). Applying Lemma [O] again we get 
£5 ■ u> = 0. Hence w = 7r(£s)n; = 0. 

The invertibility of A implies that ( Pffi ) holds with K = 3. Inequality fl2"5|) holds 
because the rank is bounded by N = 6 minus the dimension of the eigenspaces 
corresponding to eigenvalues not enclosed by 7 + which is 3. 

We can now apply Proposition |3l^. We find, microlocally near 7, a parametrix 
U for the initial boundary problem. Hence 7 ^ WF (A/ — (y ■ S(U f))\Rxdn) if 
WF / is contained in a small conic neighbourhood of 7. Here S(u) denotes the stress 
tensor which corresponds to the displacement u. The displacement-to-traction map 
it 1 — > v ■ S(u) is a first order differential operator with principal symbol 

(39) s(x, := X(u ® £) + ® + M£ • Id 3 " ^) Ids • 

It then follows from Proposition |3.3| that Ae^ 1 in a conic neigbourhood of 7. 

It remains to prove the formula ( j36[) for the principal symbol of A. From (|3^) we 
obtain 



(40) A x v 



z s A v G £5 if w G Cp, 
zpAju G C£ P if v G C£ 5 . 



From fl2~T|) we get a formula for the principal symbol of the normal derivative —D V U 
followed by restriction to M x dVL. Using this, ([39]), and ( f40"D we deduce fl36|). □ 

Remark 4. Given L > we choose S > suc/i i/iai T"*(R x dCl) \ T$ C £s holds for 
every operator in C(L, 1/2). We then choose < e < 1/2 suc/i i/jai /or every P G 
£(L,e) the initial boundary problem with Dirichlet boundary conditions is well-posed 
and the assertions in Proposition \^. 6 /\ hold. It follows from the proof of Proposition 
that microlocal forward and backward parametrices exist in the hyperbolic, in the 
mixed, and in part of the elliptic region. 

5. Propagation of Polarization 

Let L > 0. Choose < 5, e as in Remark |] at the end of section |j. In the following 
we assume P G C(L,e). 

We analyze the polarization of solutions of Pu = by applying polarization filters, 
i.e., certain approximate projection operators, to the Cauchy data of u. Fix E G 
^ fl (R x dVt) scalar elliptic with principal symbol e. The Cauchy data of a solution 
of Pu G C°°(M x Q) are, by definition, Cu = (Eul^xdn, v • S\^qq). Notice that 
the Cauchy data of a solution u of the initial boundary value problem Pu = 0, 
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u — f at R x dQ, and w = initially, are represented using the DN map as follows: 
Cf = (Ef,Af) = Cu. Here, abusing notation, we also defined Cf. 

We now describe, on the principal symbol level, the spaces onto which polarization 
filters project. Let 7 = (t, x; r, £,) G T*(R x dtt) \ 0. Set 

B$ /P (l) = ( s ( X f-4 /p u) ) **p(t,x,r,t- z% p v) c C 6 

if 7 G TCs/p- Here and are the forward and backward roots of qs/p(t, x, r, £ — 
zz/) = 0, respectively Also define the linear subspaces 

Bs,p(l) = £ ( s(}T~U ) ^v[t,x,r,C- z V ) 

where the sum ranges over the roots of qs/p(t,x,r,£ — zv) = 0. Clearly, Bs/p(j) = 
Bgi p (j) + Bg,p(y) if 7 G 1~Ls/p- The disjoint unions 

B s/p = ^ieHs/P B ^/pi.l) and B S /p = U^Gs/pBs/p(l) 
are subsets of the trivial bundles, C 6 . 

Lemma 5.1. Bf, resp. Bp, are vector subbundles ofC 6 overHs, resp. Hp, of ranks 
2, resp. 1. Bp is a vector subbundle of C 6 over Tip U Bp of rank 2. Furthermore 

C 6 = 5+ © B^ © 5+ © B P over Hp, 

(41) 

C 6 = B+®B S ®B P overH s n£ P . 

Proof. Given a point in 90 we introduce coordinates x = (x\, . . . ,x n ) = (x',x n ) 
such that fl and dQ correspond to x n > and x n = 0, respectively. Let £ = (£',£ n ) 
denote the dual variables. We also arrange that, at the given point in the coordinates, 
formula (|39|) still holds. At x n = 0, 

eld 3 \ / /iild 3 \ / (r 2 + e ,2 ) 1/2 Id 3 

X, fi > and (0) imply the ellipticity of ds/d^ n . Hence 

h Id 3 
* <9s/<9£ n 

defines an elliptic 6x6 symbol of order 0. Therefore, it suffices to prove the Lemma 
when, in the definitions of the vector spaces 5^(7), the symbols e and s(x,£) are 

replaced by (t 2 +^ /2 ) 1 / 2 H3 and £ n H3, respectively. Having made this replacement the 
assertions follow from spectral decomposition of the first order symbol g associated 
with p in Lemma |3.2| and from together with the known dimensions of kerp. □ 

Let 7To/ p and t^s/p denote the projectors associated with the decompositions (pHp. 
In the following Tl s , p and Hs/p denote 6x6 systems of pseudo-differential opera- 
tors of order having principal symbols K s , p and n s /p. The operators are defined 
microlocally where the symbols are. 

We now state how polarization in solutions can be tested on the Cauchy data. 



(42) h(t,x',T,e) 
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Proposition 5.2. Letu E D'(R x fi) 3 such that Pu E C°°(R x ttf , Cu E if (s_1) (R x 
Let 7 G T*(M x dfi) \ 0. If^EH P then 7 e WF (s) (n$Cu) if and only if WF^J" 15 (it) 
contains a Hamilton orbit above the forward /backward compressional wave bicharac- 
teristic which issues from 7 mto £/ie interior. If^E Tis \ Qp then 7 E WF^(II^Ctt) 
if and only z/WF^^w) contains a Hamilton orbit above the forward /backward shear 
wave bicharacteristic which issues from 7 into the interior. 



Proof. Introduce coordinates the proof of Lemma IO 

Abbreviate D x = (D',D n ). The simplified Cauchy data 

C u= ((D t , D')u\ Xn=0 , D n u\ Xn=0 ) 

are related to the Cauchy data Cu as follows: Cu = HC u where H E \P° (M. x <9fi) 6x6 
with principal symbol equal to h of equation It now suffices to prove the 

Proposition with Cu replaced by Cqu and the operators n^ p replaced by H~ l YY^j p H . 
The assertions now follow from Proposition |2.3| if we recall the argument in Remark |2| 
of section |3|. □ 

Curves which are bicharacteristics over the interior and reflected at non-glancing 
boundary points, with or without conversion between shear and compressional mode, 
are called broken bicharateristics. The propagation of singularities in the Cauchy data 
is stated recursively as follows. 

Proposition 5.3. Let f E H ( c s) (R x dQ), s > 3, with WF (s+1) (/) cT 5 . Let T E R 
such that no forward broken bicharacteristic which issues from 

wf (s+1) (/) n (WF (s) (n+c7/) u WF (s) (n+c/)) 

intersects Q D {t < T}. Then, after intersection with {t < T}, 

WFW(C/) = WF(* +1 )(/) U S s (WF^(U^Cf)) 

1 u5 P (WF (s) (n+c/)). 

Proof. Let u denote the solution of Pu = with Dirichlet boundary value / 
and zero initial data. By |pcn82| , Theorem 4.2] and [|Cer85|| the polarization set 
WFpoT 1 ^) n {t < T} is contained in the union of Hamilton orbits which lie above 
the broken bicharacteristics which issue from WF^ 1 ^/). We apply Proposition |5\2 
at both endpoints of bicharacteristics which connect boundary points. We obtain 
5 s/P (WF (s) (n+ p C/)) C WF (s) (IT- p C7). Therefore the right hand side of © is 

contained in the left hand side. To prove the opposite inclusion let 7 E WF^(CJ) \ 
WF (s+1) (/). Then 7 is the endpoint of a forward bicharacteristic contained in 
WF (s+1) («) and, by Proposition 0, issued from WF (s) (n+C/) or WF (s) (n+C/). 
□ 

Proposition |5.2| combined with the following result permits us to specify, without 
having to know the coefficients of P, sources for which compressional singularities 
are muted. 
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Proposition 5.4. Choose M G x dQ) 3x3 such that its principal symbol in 

equals at every (t,x,r,£\) G T*{R x dQ) the orthogonal projector onto the one- 
dimensional subspace ofM. 3 which is orthogonal to £ and v(x). Let 7 G Tip. There is 
a conic neighbourhood F GTip 0/7 such that the following inclusion holds for every 
f G H^ s) (R x dn) 3 , s>3, with WF(/) C T: 

(44) WF (s) (n+C/) C WF (s+1) (/ - Mf). 
Proof. First we show 

(45) (™ °)C«C B J + B S -. 

To see this let 7 = (t, x, r, £|) and a G C 3 with ^ ■ a — u ■ a — 0, v — v[x). In view 
of fl32l ) a belongs to the kernel of p(t, x, r, £5). In view of the definition (|39|) we have 
s(x,£f)a = {i%,v) s a- Hence (e(7)a, (ff,i/) s a) G £f (7). (f£ -£5,^)5 ^ because 
z s 7^ z s • Therefore we obtain (0, a), (a, 0) G B$ (7) + -B5 (7) proving (|45|). 

7rp = 7Tp + 7Tp vanishes on i? J + B s . Therefore ([45]) and the symbol calculus imply 

(46) n P (-;) 6 ^. 

Choose the conic neighbourhood T of 7 in such a way that the DN map A is 
a pseudo-differential operator in V x V . Shrinking V if necessary we may assume 
that every solution of Pu = which has zero initial data and Dirichlet data / 
with WF(/) C T does not contain backward bicharacteristics issuing from F in its 
wavefront set. 



Observe from formula (36) that the principal symbol of A maps the space onto 



which m projects into itself. Hence AM — MAM G \1/ and therefore 

tAn\ f EM \ (M 0\/£Id 3 \ , T0 

( 47 ) a m = n m a J ) mod *°- 



AM j ~ \ M j \ AM 



Let / G m s >(R x afi) 3 , s > 3, with WF(/) C T. Equations © and (|7|) imply 

n P ( ^ ) g fr«(R x an). 

Assume 7 £ WF (s+1) (/ - Mf). Then, recalling C/ = (Ef, Af), we obtain 
(48) 7 £ WF (s) (n P C7). 



Proposition |5.2| and our choice of F imply that ( [4"8"D holds with Hp replaced by n p . 



Hence (|48|) also holds with Hp replaced by Hp. □ 



6. Proof of Theorem 1.1 



The idea is to recover the lens maps from the elements of WF(A/) \ WF(/) with 
least time where / ranges over point sources. 

Let L > 0. Choose < 5, e as in Remark f| at the end of section £|. Let P (1) , P (2) G 
C(L,e). Assume A^) = A^ 2 \ We show that the shear and compressional lens maps 
are equal: S§ = Sf^ and Sp^ = Sp\ Let S denote the union of the sets SgJ P fl 
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((Q (k) x T 8 ) U (T s x g^)) . Observe that S^J P \S is open and dense in S$ p . We first 

prove that the lens maps agree outside S. 

Fix s > 3. Below we choose, given 7" 1 G Ts nonglancing, sources / with the 
properties 

(49) / G H^ S \R x dQ) 3 , WF (s+1) (/) = M +7 in . 

We study the singularities of the Cauchy data Cf = (Ef,A^f) = (Ef,A^f). 
Because of the zero initial condition the backward bicharacteristics issuing from 7 111 
are disjoint from WT (s+1 \u). By Proposition |J 7 in ^ WT (s) (U~ ( /P ] C f). 7 in G 

WF (s) (C/) by assumption on / since E is elliptic. Therefore 7 in G WF (s) (nJ (i) C/) U 

WFW(n+ (j) c/)u4 3) - 

Let ( 7 out ,7 in ) g S [ s ] \ S. Choose / with (g|) and 7 in £ WF (s+1) (/ - Mf). 
Proposition [5]4] implies 7 111 ^ WF^(np" (7/) for j = 1 and j = 2. Consequently, 
7 in G WF (s) (n+ (i) C7). Proposition |1| with pW implies 

(50) WF (s) (C7) n {t < t(7 out )} = M+7 in U M+7 out . 

(|50f) and Proposition |5.3| with p( 2 ) imply (a7° ut ,7 in ) G 4 for some a > 0. The 
covariable r is constant along bicharacteristics. Therefore a = 1. Thus we have 
shown 4 \<S C 4 ■ Interchanging P* 1 ** with P^ we obtain 

(51) 4 1) \5 = 4 2) \5. 

Let (7 out , 7 in ) G 4^ \ S. Choose / with (g|) and 

(52) 7 in G WF (s) (n+ (1) C7) \ WF (s) (nJ (1) C/). 

Proposition |5l| with P^ implies (|50"D. (|50|) and Proposition |5.3| with P^ imply 
(a7 out ,7 in ) G S]P U S ( p ] with a > 0. Again a = 1 follows. If ( 7 out ,7 in ) g 5^ 2) 
then also ( 7 out ,7 in ) g 4° by (0). Hence ( 7 out ,7 in ) g S ( p } if ( 7 out , 7 in ) £ S^- 
Suppose (7° ut ,7 m ) G 4 ■ Choose / as above but now with projectors ITgjp in (|52|) 
replaced by n^fp • Equation (p0|) still follows from Proposition [5]3] with P^ 1 ) because 
( 7 out ;7 in) G n 5(1) _ yn ^ WF (s) (n+ (2) C/), (0), and Proposition |^ with P( 2 ) 
imply (7° ut ,7 in ) G <Sp . Thus we have shown S P \<S C «Sp . Interchanging P^ with 
P^ 2 -* we obtain 

(53) 4 1J \ $ = 4 2) \ S- 

It remains to show that (|5l|) and (|53| ) hold with 5 replaced by the empty set. 
Assume (7° ut ,7 m ) g 4 1 . We use a limit argument to prove (7° ut ,7 in ) g 4 . First 
we observe that 7 in ,7° ut ^ £/£\ Suppose not. Then a neighbourhood of (7° ut ,7 m ) in 
4 contains a point which is in (Ts x 4^) \ S or in (S^ x Ts) \ S. This point cannot 
be in S$ . This contradicts ([51~D. Choose a sequence (7fc Ut ,7™) G 4 fl 4 ,2 ' ) > ^ ^ N, 
which converges to (7° ut ,7 in ). Let 7^ denote the shear wave bicharacteristic for P^ 
with 7fc(0) = 7™ and jk(tk) — 7fc Ut , h > 0. The length sequence of the corresponding 



18 



SONKE HANSEN AND GUNTHER UHLMANN 



sequence of geodesies is bounded. By compactness there is a limit geodesic and thus 
a bicharacteristic 7 : [0,T] -> T*(R x H) of with 7(0) = 7 in and j(T) = 7° ut . 
It suffices to show j(t) lies over the interior when < t < T. Suppose we had 7(f) 
above the boundary for some < t* < T. Then every neighbourhood of (7(t*),7 m ) 
has non-empty intersection with hence, by (51), also with Sg . This contradicts 



the continuity of the map at 7 111 . Hence we have shown Sg C . The other 
inclusions are proved in the same way. 

7. Proof of Lemma [O] 

Let 7 = (t, x, t, £1) G r<5 \ Q. To ease notation we drop the coordinates (£, x). 
Recall the definitions of the symbols qs/p(T,£)/p = r 2 — |£||/ P , the metrics p\£\ P = 

P\€\s + ( A + ^)l£| 2 > = ^l£| 2 + R£ • €, and the characteristics £ s /p = £ - ^s/p^, 
qs/p(T,£ s /p) = 0. We have |z/| = 1. 

The equation qsij, £ — zv) — qp(r, £ — zv) = (A + — zp) 2 holds for z G C. It 
implies £s • £p = £|/p = if zs = Zp. Also it implies qpij, £5) = if £| = 0. So if we 
had £f = then z s £ K because of the real principal type property of gsqp. Since 
Im zs/p > this can only happen if zp = zs- Therefore £f = implies £s ■ £p = 0. In 
the same way £p = implies £s • £p = 0. Therefore, it suffices to prove the inequality 

m- 

Choose < e < 1/2. e > will be decreased further depending on L and 5 only. 
The smallness assumption (||]) on the residual stress tensor implies 

(54) l-e<^t, , 2 <1 + ^ if^O. 

ml 2 (A + 2//)|?7| 2 

Consider the elliptic case, 7 G £5 C £ P . Without loss of generality we assume 
£ • 1/ = 0. Then £ 5 • £ P = |£| 2 + z S 2 P and |£| = |£,| < \r\/5. (A + 2//)/p < L 2 by (§. 
Therefore, we assume 

(55) pr 2 >2£(A + 2/i)|£| 2 >2^|£| 2 . 

z = zs/p is the solution with positive imaginary part of the quadratic equation 

\v\ 2 s/P z 2 - 2bz + (|£|| /p - r 2 ) = where b = R£ ■ v. 

Notice that the signs of the real parts of zs and zp are equal to the sign of b. Hence 
zszp £ R and thus £5 • £p 7^ if b 7^ 0. Assume 6 = 0. We solve the quadratic 
equations and then estimate using ([54]) and (|55|): 

|2 _P(l£l!-r 2 ) Pd^l^-r 2 ) 



|^P| 



< 



pIHI pW\p 

(p(l + e) - 2ep) |£| 2 ((A + 2p)(l + e) - 2s(\ + 2p)) |£|< 



(l-e)p\v\ 2 (l-e)(X + 2p)\v\ 2 

< lei 4 - 

Hence |£| 2 + z s z P ^ 0, i.e., (0) holds. 

In the mixed case, 7 G £p r\7is, we have G K and Imzp > 0. This implies (ff7|). 
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Consider the hyperbolic case, 7 G Tip C Tis- Without loss of generality we assume 

(56) (M P = 0. 

Then the roots of = r 2 — |£ — zu\ P = — \u\pz 2 — |£| p + r 2 have opposite signs. 
Since qp < qs this is also true for the roots of = r 2 — |£ — zv\ 2 s . Furthermore 
< \zp\ < \zs\- Since zs and zp are both forward they have the same sign. For 
simplicity we assume < zp < Zs- We shall use the estimate 

(57) £s-£p> \i\ 2 -ez s \i\ + z s z P . 

(A + 2/i)£ • v + • v = is equation (|56|) restated. From this and @ we deduce 
2|f • v\ < e|£|. Hence |(z s + 2p)(f • f)| < 2;z 5 |£ • i/| < ez s \£\. (0) follows. 
The equation 95 (r, £ — 2:5^) — qp(r, £ — ;zpp) = is equivalent to 

(58) - *>|p) - 2z 5 (£, - ((A + ^)/p)\i\ 2 = where t = z P /z s . 

We estimate the root z s of this quadratic equation. Using the Cauchy-Schwarz in- 
equality and fl54|) to estimate \u\s from below we deduce from 



N/4< |f|s + |flV(A + //)/// H2t\v\ p < 



v\ 



s- 



Decreasing e > if necessary, we assume ezs < |£| if 2£|z/| p < \v\ s . Inserting 
this estimate into ( [ST]) we get £s • £p > z s z P > 0. It remains to prove ( j37l) when 
2t|p| p > \u\ s . From ([54] ) and (f|) we get |^|p/|^| s < 3(A + 2/x)//x < 3L 2 . Decreasing 
e > if necessary, we assume e 2 \v\ p < 2|p| 5 . Hence e 2 z s /A < zsZp. We estimate the 
right hand side of fl57f) from below and get £s • £p > (|£| — ^s/2) 2 > 0. 
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